Numerical solution of Laplace's equation